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Much of the energy of the nuclei colliding at RHIC or SPS is converted into final state 
^ I hadronic particles. About a quarter of this energy is in baryons and antibaryons. There 
are nearly 10 strange quark pairs per central rapidity participant. Do we really understand 
O ■ the hadronic particle yields? Do we need to introduce post-Fermi-model ideas such as 
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chemical non-equilibrium in order to understand how a deconfined state hadronizes? 
1. STATISTICAL HADRONIZATION 



Statistical Fermi- Pomeranchuk models have been used extensively to study particle 
O I yields and spectra since 1950 [|1],0. This approach was developed as a qualitative de- 
^ I scription of the gross features of particle production. It was originally not meant to be 
theoretically an accurate picture, just a phase space estimate of what Fermi called an 



upper limit on particle production in strong interaction processes. In fact, Fermi defined 
the limitations of the statistical approach stating three conditions: 1) exclusion of par- 
', tides such as photons which have weak coupling to the interacting system; 2) absence 
of absolute chemical (abundance) equilibrium for many semi-weakly coupled particles; 3) 
requirement of relative baryochemical equilibrium. At that time strangeness was not yet 
discovered, and today we can add 3b) requirement of relative strangeness equilibrium. 

In the ensuing decade, another important feature of the strong interactions was discov- 
ered: the existence of numerous hadronic resonances. This property of hadronic interac- 
tions poses a challenge for the statistical hadronization model as the yield of particles is 
sensitive to the unidentified high mass hadron resonance states. We will discuss how this 
influences the expected particle yields in section |^. Our objective is to establish the mag- 
nitude of systematic theoretical error we have to expect given incomplete knowledge of the 
hadronic mass spectrum. A tacit assumption in our approach is that hadron-hadron in- 
teraction is well described by the formation of resonant states, and the remaining residual 
forces are negligible in comparison to the available energy content per particle. 

Understanding the hadro-chemistry, i.e., of the composition of particles produced, can 
flnger-print the phase of matter which has undergone statistical hadronization. Therefore, 
we will describe how chemical parameters, the chemical potentials fii and the phase space 
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occupancy parameters 7^, allow to regulate the particle number. We show how the chemi- 
cal non-equilibrium is characterized in section |^. We then demonstrate, in section ^ that 
the nearly complete particle production results force upon us the study of chemical non- 
equilibrium. We demonstrate in particular that a precise description of all RHIC results 
obtained at a/snn = 130 GeV (RHlC-130) is possible within the scheme we propose. 



2. EXPONENTIAL HADRON MASS SPECTRUM 

Fermi's statistical model provided the intellectual environment for the theoretical un- 
derstanding of the exponentially growing hadronic mass spectrum within the statistical 
bootstrap model invented by Rolf Hagedorn [Q. Hagedorn recognized early on that the 
large number of different hadronic states, and the increase in their number with their 
mass, has great practical implications for the behavior of matter at high temperature, 
and can lead to the formation of a new phase of matter. 

Within Fermi's statistical model the number of heavy hadron resonances produced is 
exponentially suppressed with the ratio of the mass of the particle considered to the 
temperature present. In spite of this, if the number of different resonances per unit of 
mass, the mass spectrum p{m), increases exponentially. 
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there still would be a critical temperature of hadronic matter, corresponding to the (in- 
verse) slope Th of the exponential mass spectrum: 
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Hagedorn called this the boiling of hadronic matter . Hot hadronic matter behaved at 
this point in the same way as in high temperature 1st order phase transition |^,^. 

There is strong empirical evidence that the experimental mass spectrum grows exponen- 
tially at least in the interval of interest to the study of particle production . Smoothed 
mass spectrum based on 4627 states known in 1996 is shown as a short dashed line in 
figure while the 1411 states of 1966, used by Hagedorn 0, are depicted by long dashed 
line in figure |I]. The solid line is a fit to the experimental mass spectrum Eq. (|I|) with the 
assumed value a = —3 which yields Th = 158 MeV. 

Proposing that the solid line describes reasonably the mass spectrum, we see that a 
significant number of resonances remains to be identified for m > 1.4 GeV (note the 
logarithmic scale in the figure |l|). When a study of particle yields is made, a significant 
systematic uncertainty derived from the lack of knowledge of hadronic mass spectrum can 
arise, in particular so for models finding relatively high statistical hadronization tempera- 
tures. At a hadronization temperature of Th = 145 MeV, a study of the properties of the 
particle yields suggests that only about 5% of pions are 'missing'. However, at Th — 175 
MeV []TU[, this increases to 35%, and the high value of temperature considered poses a 
practical challenge, since the Statistical Bootstrap model and the Lattice-QCD results for 



2-1-1 flavors [|TT| are converging to a critical temperature at Th < 160 MeV. 

One way to assess the magnitude of the systematic error is to compare the properties 
of the hadron gas using the mass spectrum of the type given in Eq. (P with that studied 
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using the currently known experimental mass spectrum. For a range of values a = —2.5 
(most divergent curve in left frame in figure ^ and a = —7, the energy density e (solid 
lines) and pressure P (dashed lines), both scaled with T^, are compared to the values 
obtained summing all known hadronic states (thin line in figure Q). The differences as 
expected are dramatic at high temperature, near to the boiling point of hadronic matter. 

It turns out that the problem is reduced considering the hadron proper volume 0. 
Reconsidering the energy density assuming a proper energy density of hadrons as derived 
from the bag constant i3=(190 MeV)^, we find that the expected excess counting the 
energy density degrees of freedom remains at the level of 10% as is seen in right frame in 
figure 0. While this allows to apply statistical hadronization method in study of particle 
abundance ratios, the influence of quantitatively uncertain finite size correction is very 
large, implying in particular in the pion yield a considerable systematic theoretical error. 

3. CHEMICAL NONEQUILIBRIUM 

In general the fugacity Tj of each individual particle will comprise the two chemical 
factors associated with the two different chemical equilibria. For example, let us look at 
the nucleon, and the antinucleon fugacity and chemical potentials: 

Ttv = 77ve^^/^, = 7;ve-'^B/T. = + Tln7;v, aj^ = -fiB + Tln^N.{3) 

There is an obvious difference between the two chemical factors in Eq. (0): the number 
of nucleon-antinucleon pairs is associated with the value of 'y^ but not withyUs- This can 
be seen looking at the first law of thermodynamics, in this context written as: 

dE = -PdV + TdS + aNdN + ajjdN, 

= -PdV + TdS + fiBidN-dN)+T\n-fN{dN + dN). (4) 

To obtain the second form, we have employed Eq. (^. We see that /ib is the energy 
required to change the baryon number, B = N — N, hj one unit, while the number of 
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Figure 2. The energy density e (solid lines) and pressure P (dashed lines), both scaled with 

T^. Left: point hadrons with an exponential mass spectrum (thick lines) arc compared to results 
using known hadrons (thin line). Right: hadron gas with finite volume correction, see text for 
more detail. 



nucleon-antinucleon pairs, 2iVpair = iV + A^, is related to 7Ar. For 7jv — ^ 1, the last term 
vanishes, at this point small fluctuation in number of nucleon pairs does not influence the 
energy of the system, we have reached the absolute baryochemical equilibrium. 

It is convenient to follow the quark flavor even in the study of hadron yields, since this 
allows to keep the same notation across the phase boundary of quark matter and hadronic 
gas matter. We use: 

A.^e'^"/^, Xa^e^^/^, A, = e'^^/^ = + ;,,), A^ = A„A,. (5) 

When the light flavors -u, d remain indistinguishable we introduce Ag. 

The relationship of baryochemical potential /xb and of strangeness chemical /xs to quark 
chemical potential is: 

^x■B, = 3//g, Ab = Aj, As = ^, //s = J/^B - /^s- (6) 

As 6 

These relations arise considering that three quarks make a baryon, and remembering that 
strange quarks carry negative strangeness and one third of baryon number. 

When using of fugacities which follow the valance quarks, Aj, i — u,d, s, hadronic 
particle yields can be easily checked, and thus errors and omissions in a rather complex 
and large particle 'zoo' minimized. Consider, as an example, the ratio E~ (dss) /E~ (dss) . 
Given the quark content and ignoring (at RHIC very small) isospin asymmetry, we find: 

^ XiXq 2 4^ 

Since all S resonances which contribute to this ratio are symmetric for particles and 
antiparticles, and possible weak interaction feed from fl{sss), and respectively fl{sss), 
are small, these expressions are actually nearly exact. Thus we have inverted the relation 
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and expressed /is in terms of hb and the cascade ratio. This allows the reader to check 
for correctness results presented elsewhere. Results shown in [|1^] appear inconsistent by 



50% while those in |12| are numerically consistent. 

The important lesson to be drawn from Eq. (|^ is that when we compare in a ratio 
particle with antiparticle we 'see' the chemical potentials. As we shall see below and use 
in section ^, the chemical non-equilibrium parameters are probed with products of yields 
of particles and antiparticles. 

Near to chemical equilibrium, we use three non-equilibrium parameters, 7^, and 7^, 7^ 
or equivalently just two parameters introducing, jq = y/juld- In quark matter, these 
coefficients express the approach to the expected chemical equilibrium yield by the quark 
abundances. In a coalescence hadronization process, gluons fragment into quark pairs and 
the net yields of quarks and antiquarks of all flavors is redistributed among all individual 
hadrons. Hence even if there were no change of the quark pair number in hadroniza- 
tion, the values of 7^, 7^, 7^ in hadron gas and quark matter must differ. Moreover, the 
phase spaces have different size, and it is impossible in the rapidly evolving fireballs to 
reequilibrate several quark fiavors. Thus, we have to distinguish: 

7.^^^7.^^^7?^^ from 7^,7.^^,7^- 

The lesson is that in an analysis of experimental data, we explore the properties and 
parameters of the hadron gas phase (HG). While these may be smooth across hadroniza- 
tion for the chemical potentials fii, we expect the phase space occupancy parameters 7^ 
to be quite different in the confined and deconfined phase. We obtain the fugacities of all 
hadrons in terms of six parameters, which implicitly carry the upper index 'HG'. 

It is best to see how this works looking at typical examples: 

a) baryons: protons p{uud) , antiprotons p('U'U(i), A{uds), Q{sss), etc: 

Tp = 7'7de2'^-+'^^ Tp- = 7^7de-2'^"-^^ Ta = 7«7d7. e'^"+'^''+'^% = 7^ e-3^^ (8) 

b) mesons: tt^ (ud), -K^iud), K^{us), 4>{ss), etc: 

= 7«7de^--^^ T^- = 7«7de-^"+^^ T^- = 7.% 6"'^"+'^% = ^1 (9) 

Note that in the products of particle Tp and antiparticle Ta fugacity the chemical po- 
tentials fii always cancel TpTa = /(7j)- 

In the statistical hadronization approach, we have the same value of T for all hadrons 
with the same valance quark content. For example Tp = T^+. Thus, we assume tacitly 
that the relative population of heavier resonances is in chemical equilibrium with the 
lighter states. This means that statistical hadronization resolves the valance quark and 
antiquark distribution and does not allow for the possibility that heavier resonances may 
simply not be populated. This method, thus, is most suitable for a hadronizing quark 
matter fireball, and may miss important features of a hadron fireball which never entered 
the deconfined phase. 

The physical properties of the gas of hadrons with 7„, 7^, 7^ 7^ 1 have certain highly 
desirable properties. Consider the properties of the pion gas as function of 7^, and in 
particular the entropy (B for Bosons and F for fermions): 

^^/^= = , i , • (10) 

J [zirn) 7^ e r — 1 




Here, E^^ = \/m?^ + jp-. We see, in figure |], that the entropy density is rising rapidly and 
it nearly doubles with 7^ increasing from the equilibrium value 7^ = 1 toward condensa- 
tion point 7q = e^^l'^ . This high entropy content of the oversaturated pion gas allows 
the sudden hadronization of QGP without reheating or infiation (volume expansion). A 
striking feature of the experimental data analysis which allows for 7^ is the maximization 
of entropy density in pion gas, i.e., natural tendency toward the value 7^ = e^^l'^ . 

4. DATA ANALYSIS 

There are in our approach at most 6 chemical parameters, and hadronization temper- 
ature. Especially at RHIC, we need not distinguish the light fiavors m, d and thus the 
number of statistical parameters is reduced by two. A further reduction is arrived at 
when we demand that strangeness balances anti-strangeness locally, and finally the re- 
quirement that when entropy rich quark-gluon plasma hadronizes, the yields of hadrons 
are maximizing the entropy content of hadrons, yields 7^ = e^^l'^ . Thus, we are in fact 
needing just two chemical parameters which are usually \q and 7^. Both at SPS and 



at RHIC |]T^, there is good evidence for a single freeze-out to apply, in consistency with 



the sudden hadronization hypothesis. In other words, one single temperature allows to 
understand both yields and spectra of hadrons. 

The computation of the particle yields is much simpler than for SPS. For central ra- 
pidity, we have, at RHIC-130, approximate longitudinal scaling. Thus, we can act as 
if a series of fireballs at all rapidities was present. Then, we simply evaluate the full 
phase-space yields in order to obtain particle ratios. We do not include in our analysis 
trivial results such as 7r"'"/7r~ = 1. We also do not fit the results for K* and K* since the 
reconstructed yields depend on the degree of rescattering of resonance decay products. 

We first convince ourselves that the introduction of chemical non equilibrium in the 
study of hadron abundances is necessary — even though one could argue that if it is not 
needed, the data fit will converge to chemical equilibrium condition. Therefore before we 
present an overall fit of many particle yields, we present a short but persuasive argument. 

As noted, when we evaluate the product of particle and antiparticle yields, the chemical 
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Figure 4. Ratios of products of particle and antiparticle abundances in the T-^sMq plane. 
Ranges derived from data shown in table [l|. 



potentials Hi cancel. Studying such products allows to focus only on T, 7^,7^. We next 
identify three ratios involving products of particles and antiparticles: 
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By comparing mesons with mesons, and baryons with baryons we reduce uncertainties 
about excluded volume. 

The ratios considered in Eq. (^) have only two parameters T and 7^/75. Thus, it is 
possible to present, in a two dimensional figure, how these three ratios behave, as is seen in 
figure ^ for the experimental data obtained at RHIC-130 (see table |I] below). Dashed area 
shows the allowed parameter domain. The cross in the figure is the result of global data 
analysis discusses below where the ratio of baryons to mesons, not used here, fixes the 
temperature. While the kaon to pion ratio would tolerate within one standard deviation 
the chemical equilibrium, the baryon double ratios, and in particular the more strange 
cascades, are clearly demanding a value of 75/7^ > 1. 

We use the latest experimental results for the ratio of hyperons to nucleons (yields of 
p, p are weak decay feed corrected. A, A are uncorrected) [|16|. This result is nearly by 
a factor two different from the experimental data stated in ||12[, and this reference also 
does not consider S, S results other than their ratio, and is therefore consistent with the 
chemical equilibrium in the result of its analysis. 

A global fit to the experimental data for RHIC-130 is given in table |1[ We consider 
here 21 particle ratios. In some cases within the error other results from same and/or 
another collaboration are available, our selection is somewhat subjective, but does not 
influence in essential way the general conclusions which follow. In the three last columns, 
the results for both chemical equilibrium (last column) and non-equilibrium fits are seen. 

Next to the fitted ratios, we show in parenthesis the contribution to the error (x^) for 
each entry. We consider statistical errors for the experimental results, since much of the 
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Table 1 

Central-rapidity hadron ratios at RHIC-^snn = 130 GeV. From top to bottom: experimental 
results, fitted chemical parameters, the physical properties of final state hadron phase-space, 
and the fitting error. Columns: ratio considered, data value with reference, two non-equilibrium 
fits, and in the last column, the chemical equilibrium fit. The superscript * indicates quantities 
fixed and not fitted. The superscript ^ indicates the error given is dominated by theoretical 
considerations. Subscripts H, A mean that these values include weak cascading given in heading 
of table. In parenthesis, we show the contribution of the particular result to the total ^ . 



Data 



Ref. 



100% H ^ y 

40% Y 



40% H ^ y 
40% Y 



40% H ^ y 

40% y ^ iv 



p/p 

As /As 

K-/K+ 

K-/7r± 

K+/^± 

As/h- 

As//i- 

As/p 

As/p 

E-fTT 

E-/h- 
E-/A 

n/E- 
H/eF 

^}/h- 

n/n 

<p/K- 

<p/h- 



0.71 ± 0.06 
0.71 ± 0.04 
0.83 ± 0.08 
0.87 ± 0.07 
0.15 ± 0.02t 
0.17 ± 0.02t 
0.059 ± 0.004t 
0.042 ± 0.004t 
0.90 ± 0.12 
0.93 ± 0.19 

9.5 ± 2 
13.4 ± 2.5 
0.0088±0.0008t 
0.0085±0.0015 
0.0070±0.001 
0.193±0.009 
0.221±0.011 



0.0012±0.0005 
0.95±0.1 
0.15±0.03 
0.021±0.001 



PI 
17| 

PI 

PI 
18 



PI 
PI 



[16| 
PI 



|18| 
PI 



[17| 



0.672(0.4) 
0.750(1.0) 
0.793(0.2) 
0.925(0.6) 
0.156(0.1) 
0.169(0.0) 
0.057(0.6) 
0.042(0.0) 
0.805(0.6) 
0.899(0.0) 
9.4(0.0) 
13.7(0.1) 
0.0092(0.2) 
0.0076(0.3) 
0.0061(0.9) 
0.190(0.1) 
0.207(1.6) 
0.20 
0.22 
0.0015(0.4) 
0.87(0.7) 
0.174(0.6) 
0.022(1.3) 



0.677(0.3) 
0.747(0.8) 
0.803(0.1) 
0.924(0.6) 
0.157(0.1) 
0.170(0.0) 
0.049(6.7) 
0.036(2.0) 
0.662(3.9) 
0.731(1.1) 
9.2(0.6) 
13.3(0.0) 
0.0097(1.2) 
0.0079(0.2) 
0.0064(0.4) 
0.189(0.2) 
0.206(1.9) 
0.21 
0.23 
0.0016(0.7) 
0.88(0.5) 
0.177(0.9) 
0.023(2.5) 



0.688(0.1) 
0.757(1.4) 
0.818(0.0) 
0.933(0.8) 
0.150(0.0) 
0.161(0.2) 
0.047(9.6) 
0.035(2.8) 
0.494(11.5) 
0.543(4.1) 
7.4(27.7) 
10.6(9.6) 
0.0069(5.8) 
0.0056(3.8) 
0.0046(6.0) 
0.132(45.4) 
0.144(48.4) 
0.18 
0.20 
0.0010(0.13) 
0.90(0.3) 
0.148(0.0) 
0.018(10.2) 



^HG 
\ 

fiB [MeV] 

^HG/^HG 

A. 

[MeV] 



140.3 ± 1.1 
1.64* 
1.0700 ± 0.0076 
28.5 
1.50 ± 0.04 
1.0243* 
6.1 



142.5 ± 1.2 
1.63* 
1.0686 ± 0.0076 
28.4 
1.48 ± 0.04 
1.0218* 
6.4 



165.8 ± 2.2 

1* 

1.0654 ± 0.0082 
31.3 

1* 

1.0186* 
7.4 



E/b [GeV] 


34.7 


34.3 


34.1 


s/b 


9.5 


9.3 


7.0 


S/b 


233.4 


227.7 


238.5 


E/S [MeV] 


148.7 


150.5 


143.0 


xVdof 


10/(21 - 3) 


25/(21 - 3) 


188/(21 - 2) 



J.R./J.L. 



26.10.2002 
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systematic error should cancel in the particle ratios. However, we do not allow, when 
pion multiplicity is considered, for errors smaller than ~ 10%, which is our estimated 
error in the theoretical evaluation of the pion yield due to incomplete understanding of 
the high mass hadron resonances. Some of the experimental results are thus shown with 
a 'theoretical' error. When such an enlargement of the experimental error is introduced, 
a dagger as superscript appears in the experimental data second column in table |I|. 

The high yields of hyperons require significant (30-40%) yield corrections for unresolved 
weak decays. Some experimental results are already corrected in this fashion: the weak 
cascading corrections were applied to the most recent p and p results by the PHENIX 
collaboration [Q, and in the S/A and H/A ratio of the STAR collaboration we use here 
ITf . However, some of the results we consider are not yet corrected |T^JT^, and are 



indicated in the first column in table |I| by a subscript A or S. We present two non- 
equilibrium fits, left with 100% E ^ Y cascading acceptance and with 40% Y ^ N 
acceptance. Then, a non-equilibrium fit with 40% S — > y and 40% Y N, and in the 
last column, the chemical equilibrium fit with 40% cascading. 

Below the fit results, we show the statistical parameters which are related to each fit, 
and at the very bottom the of the fit. We observe a considerable improvement in the 
statistical significance of the results of chemical non-equilibrium fits. The results, shown 
in the table |l], are obtained minimizing in the space of 3 parameters: the chemical freeze- 
out temperature T, and 2 chemical parameters A^, 7^, the value of 7^ is set at its maximal 
value, 7^ = e™'"/^, and the value of is derived from the strangeness conservation 
constraint. Freeing these parameters does not alter the results, the fit converges to local 
strangeness neutrahty, within a few percent, and to full pion phase space saturation. 

We note that several particle yields are not properly described in the last column of 
table and hence a large results in this chemical equilibrium fit. However, on a 
logarithmic scale only results involving A, A would be clearly visible as a discrepancy. 
The second and third last column show result of chemical non-equilibrium fits bracketing 
the S-cascading, and yielding an excellent confidence level. 

In the bottom of table ^ we see that the chemical non-equilibrium fit specific strange- 
ness content s/b ^ 9.5 is nearly 40% greater than the chemical equilibrium result. This 
originates in js/lq — 1-5 {i.e., 7^ ~ 2.5). This specific yield of strangeness and strange- 
ness occupancy we measure in the phase space after hadronization is consistent with the 
expected QGP properties before hadronization with nearly saturated strangeness phase 
space. The Wroblewski ratio p2|, W = 2{ss) / {{uu) + {dd)), is nearly proportional to 
1s/lq and it increases by 40% at RHIC-130 compared to SPS. However, an analysis with 
7^/7^ ~ 1 will not report this enhancement p3|1 . 

Our chemical nonequilibrium results yield a low freeze out temperature T/j = 140- 
142 MeV. This is considerably less than the temperature of phase transformation, Th 
which might have been naively expected. This range of hadronization temperature now 
seen agrees better with the expectations we had upon consideration of the effect of the 
fast expansion of QGP We believe that the wind of expending quarks and gluons 



adds to the thermal pressure and the combined kinetic and thermal pressure of exploding 
quark-gluon plasma can press out the confining vacuum even at lower temperature p5 . 



The results we find strongly favor chemical non-equilibrium description of the hadroniza- 
tion process. A non-equilibrium chemical analysis of heavy-ion particle yields offers pro- 
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found insight into the physical properties of the dense hadronic matter formed in the 
relativistic heavy-ion colhsions and allows to infer conditions prevailing in the decon- 
fined phase. While the equilibrium Fermi model provides a first impression about the 
range of freeze-out temperatures, our post-Fermi model resolves the 'fine structure' of the 
hadronization yields. 
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